A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

06 (08.02.2, 08.02.4, 08.09.2 (SS Lac), 10.15.2 (NGC 7209)) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OO 



Spectroscopic orbit of the ex-eclipsing binary SS Lac 
in the young open cluster NGC 7209 

Lina Tomasella^ and Ulisse Munari^ 

^ Osservatorio Astrofisico del Dipartimento di Astronomia, Universita di Padova, 1-36012 Asiago (VI), Italy 
Osservatorio Astronomico di Padova, Sede di Asiago, 1-36012 Asiago (VI), Italy 

Received date ; accepted date 



Abstract. The no-longer-eclipsing system SS Lac in the 
young open cluster NGC 7209 has been recently an- 
nounced to show a constant radial velocity. Puzzled by 
this finding, we have monitored the system during 1997 
obtaining 24 Echelle-|-CCD spectra over 8 orbital revolu- 
tions. Our spectra reveal a nice orbital motion with pe- 
riodic splitting and merging of spectral lines from both 
components. 

An accurate orbit has been derived, together with in- 
dividual masses of the stars. SS Lac presents a moder- 
ately eccentric orbit and a probably full synchronization 
between stellar rotation and orbital revolution. 

Key words: Binaries: eclipsing - Binaries: spectroscopic 
- Stars: individual (SS Lac) - Open Clusters: individual 
(NGC 7209) 
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J> 1. Introduction 



According to recent photometry (cf. Vansevicius et al. 
1997) SS Lac belongs to the young open cluster NGC 
7209 . Hoffmeister (1921), Dugan & Wright (1935), Wach- 
mann (1936), Nekrasova (1938) and Kordylewski et al. 
(1961) described SS Lac as an eclipsing binary of 14.4 
days period, with significant eccentricity because the equal 
depth secondary eclipse occurred at phase 0.57. However, 
more recent observations by Zakirov & Azimov (1990), 
Lacy (1992), Mossakovskaya (1993) and Schiller et al. 
(1996) show that the eclipses no longer occur. The end of 
the eclipsing phase is set around 1940 by Mossakovskaya 
(1993) and around 1960 by Lehmann (1991). The latter 
ascribes the rotation of the plane of the orbit to the pres- 
ence of an unseen third body in the system, orbiting at a 
great distance the closer central pair. 

Quite puzzling has been the report by Schiller et al. 
(1991) and Schiller & Milone (1996) that 1983-84 spectra 
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of SS Lac did not reveal the expected double-lined na- 
ture of the spectra. Very little, if any, variation in radial 
velocity was observed. Using available predictions for the 
rate of change of the SS Lac orbital inclination, Schiller 
& Milone (1996) expected a semi-amplitude of the radial 
velocity curve of 150 km sec~^, far in excess of their instru- 
mental resolution. They suggested that SS Lac could be 
a triple system suddenly become chaotic or that a close 
encounter with another NGC 7209 member could have 
ionized the binary (cf. Schiller 1996). 

Stimulated by the Schiller et al. (1991) report, Etzel 
et al. (1996), Stefanik et al. (1996) and Etzel & Volge- 
nau (1996) announced in lAU Circulars that, according 
to their preliminary spectroscopic observations, SS Lac is 
a double-lined system, with indications of variability in 
the radial velocities. So far no investigation of the orbital 
motion of SS Lac has appeared in literature. 

As an ex-eclipsing binary in a well populated young 
open cluster, SS Lac clearly deserved further investiga- 
tions to clarify the whole issue. In this note we report 
about our spectroscopic monitoring of SS Lac performed 
to the aim of confirming the binary nature and to derive 
the spectroscopic orbit. 

2. Observations 

Twenty-four Echelle+CCD spectra of SS Lac have been 
obtained with the Asiago 1.82 m telescope during 1997. 
They cover the range from 4300 to 6600 A with a res- 
olution of A/ A A -22,500 at H/3 (from the FWHM of 
comparison spectrum Thorium lines). The spectra have 
been reduced in a standard way with the IRAF software 
package. A sample of Ha profiles recorded on our spec- 
tra is shown in Figure 1, where the periodic splitting and 
merging of the individual components is evident. 

On our spectra we measured only the radial velocities 
of the Ha, H/3 and H7 lines (the latter over two distinct 
Echelle orders). A few other metallic lines were indeed 
visible, but too weak for meaningful measurements. The 
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radial velocities reported in Table 1 come from the hy- 
drogen lines, weighted according to the S/N of the stellar 
continuum around the line (typically 6:4:3 for Ha:H/3:H7, 
with an average S/N=50-60 around Ha). 
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Fig. 1. A sample of Ha profiles of SS Lac reckoned according 
to the orbital phase. The merging and splitting of the spectral 
lines is evident. Orbital phases on the left are according to 
Eq.(2). 



3. Spectroscopic orbit 

Our observations were secured over a time span cover- 
ing less than 8 revolutions of SS Lac and therefore are 
not suited for a very accurate period determination. A 
Deeming-Fourier analysis of data in Table 1 suggests a 
period of 14.4±0.3 days, very close to the 14.416 days 
determined from eclipse timing by several Authors (e.g. 
Brancewicz & Dworak 1980). In the computation of the 
orbit, the period 14.41638 from Eq.(2) has been adopted 
(see below). 

The orbital solution is given in Table 2. With the ex- 
ception of CI and eccentricity, the formal errors for the 
other elements are <1.5%, with a weighted deviation of 



the observations from the computed orbit of ^0.8 km 
sec~^. The solution is graphically presented in Figure 2. 

It is worth to compare the photometric and spectro- 
scopic orbital solutions (cf. Savedoff 1951) through the 
combination of the eccentricity (e) and longitude of the 
periastron (ft) elements (cf. Struve 1948): 
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The photometric value is -hO.llO (Dugan & Wright 1935), 
very close to the -1-0.108 value from the spectroscopic or- 
bital solution in Table 2 (and orbital inclination from sect. 
4.3). Thus the photometric and spectroscopic orbital so- 
lutions appear in excellent agreement. 



Table 1. Heliocentric radial velocities of the SS Lac compo- 
nents. Solid line and filled circles in Figure 2 refer to star a. The 
quoted errors are probable errors. When the two components 

were too close in radial velocity for lino splitting into individual 
components, a single value is given (central column). MJDq = 
JDq - 2450000. 
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4. Discussion 

4.1. Eclipse timing 

The MDJ0=716.315 spectrum in Figure 1 presents a per- 
fect radial velocity superposition of the two spectra, thus 
it corresponds to eclipse conditions. The eclipse ephemeris 
given by Mossakovskaya (1993) 

Mini = 2415900.76+ 14.41619(±0.00013) X S 

predicts a -1-0.456 day shift for the MJD0=716.315 spec- 
trum, which is comparable with the propagation of 
the uncertainty on the period. Imposing phase coinci- 
dence between Mossakovskaya's principal minimum and 
MJD0=716.315 spectrum, this leads to the improved 
ephemeris 

Mini = 2450716.32(±0.15)-hl4.41638(±0.00010)xi;(2) 

The slightly longer orbital period nearly coincides with 
those given by Brancewicz & Dworak (1980, P=14.4163 
days) and by Dugan & Wright (1935, P=14.41629 days). 



Table 2. Orbital elements for SS Lac. Where appropriate, the 
second value correspond to the component b represented by a 
dotted line and open circles in Figure 2. The quoted errors are 

the formal errors of the orbital solution. The entry '^deviation" 
is the mean weighted deviation of the observed radial veloc- 
ities from the computed orbital solution (weight = err~^ in 
Table 1). MJDq = JD© - 2450000. The error on the masses 
is ~0.1 M0. 



as suggested by the relative depths of the central Doppler 

cores and contribution to the overall Ha line wings in Fig- 
ure 1. A lower brightness agrees with the slightly smaller 
mass of component b (cf. the semi-amplitudes and semi- 
major axes in Table 2) and the just slightly reduced ampli- 
tude of b eclipses compared to a eclipses (a few hundredth 
of a magnitude, cf. Dugan & Wright 1935, Mossakovskaya 
1993). 
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4-2. Spectral classification 

The absence of He I absorptions, the paucity of metallic 
lines and the very strong Stark broadening of the wings of 
the Balmer lines suggest a spectral classification around 
A2 V, virtually identical for the two members (cf. Fig- 
ure 1). A more refined classification however needs de- 
voted observations at the head of the Balmer series. Com- 
ponent b appears slightly less luminous than component a, 



orbital phase 



Fig. 2. Orbital solution for SS Lac. The curves give the or- 
bital solution of Table 2. Solid curve and filled circles refer to 
component a in Tables 1 and 2. The crosses mark spectra for 
which the too small velocity separation caused merging of the 
two components into a single line. The dcished line is the bari- 
centric velocity (-21.2 km sec~^). Orbital phase according to 
Eq.(2). 



4-3. Orbital inclination 

The spectral classification of the two components can be 
turned into a physical radius of R=2.25 Rq from interpo- 
lation of tabular values given by Allen (1973). The orbital 
separation from Table 2 is asini = 0.200±0.002 AU = 
42.5±0.5 Rq. Lehmann (1991) has stated that the system 
presented total eclipses (central eclipses given the nearly 
identical size of the components) in 1900-1915, and af- 
terward the eclipses monotonically reduced in depth until 
they disappeared around 1960. From separation and size 
of the two components it may be computed that at the 
end of the eclipsing season the inclination was i ^ 83°. 
The corresponding linear change in orbital inclination is 



^=0.13 ±0.01 degyr- 
dt 



(3) 
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which is significaiitly lower than the di/dt = 0.18 deg 

yr""'^ obtained by Lchniann (1991) on the assumption of a 
B7 V classification for both components (clearly ruled out 
by our spectra that show no evidence for Hel absorption 
lines). From Eq.(3) and the occurrence of central eclipses 
at the beginning of this century, the current orbital incli- 
nation is 



*(1998) 



78° 



(4) 



The above results indicate that the non-eclipsing season 
of SS Lac lasts for about 1275±90 years. 

4.4- Rotational velocities 

The rotational velocity of the components of SS Lac has 
been derived by comparison with standard stars from Slet- 
teback (1975), which were observed under the same instru- 
mental conditions. The velocity turned out to be the same 
for both components: 



10 ± 5 



km sec 



(5) 



The dependence on inclination (see Eq[4]) has been re- 
moved under the assumption that the rotational and or- 
bital axes are parallel. A Vrot = 10 ± 5 km sec"-^ means 
that rotational and orbital motion are synchronized or 
close to, the orbital velocity amounting to Vorb = 8 km 
sec~^. 

It is worth noting that SS Lac shows evidence for syn- 
chronization when the orbital circularization is still far 
from being achieved. The two phenomena seems therefore 
to evolve over quite different time scales, in agreement 
with the theoretical expectations (cf. Zahn 1977). 

The expected synchronization time for a binary with 
the SS Lac parameters may be estimated from the Tas- 
soul's (1987) formalism to be: 



tsyniyrs) > 2500 X [P{days)]- 
> 4 million years 



(6) 



which is relatively short compared to the NGC 7209 clus- 
ter age (300 million yrs according to Lynga 1985). The 

synchronization status cannot therefore be used to decide 
if SS Lac is or not a primordial binary of NGC 7209. 

4-5. Individual masses 

A straightforward application of the Kepler's III law to 
Table 2 data gives 2.69 and 2.80 Mq for the masses of 
the two SS Lac components. The tabular mass from Allen 
(1973) for the A2 V spectral type is 2.78 M©, in excellent 
agreement. 



4-6. Mass- Luminosity relation 

Vansevicius et al. (1997) have estimated a 1.0 kpc distance 
and a Ay=0.54 mag extinction to NGC 7209. Using a 



bolometric correction B.C.= -0.29 mag for the spectral 
type A2 V (Allen 1973), the bolometric magnitude of the 

SS Lac components is 



Mboi = -f-0.25 mag 



(7) 



corresponding to a luminosity L = 65 Lq. The mean 
mass for the two components is 2.745 Mq. At this mass 
the Mass-Luminosity relation has the numerical form (cf. 
Bowers & Deeming 1984): 



L M 
log— = 0.479 -\- 2.91 log—— 

Lq ■ Mq 



(8) 



which predicts L ~ 57 Lq for the SS Lac components. The 
agreement between the observed and predicted luminosity 
is satisfactory in view of the uncertainties in the distance, 
reddening and the approximate bolometric correction. 



Table 3. Heliocentric radial velocity of the weak component 
of the Balmer linos described in Sect. 4. 7. The errors may be 
quantified as ±1 km scc^^ when N=3, ±2 for N=2 and ±4 for 
N=l. MJDq = JDq - 2450000. N = number of Balmer lines 
used. 
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4- 7. The unseen third body 

Lehmann (1991) ascribes the rotation of the plane of the 
orbit to the presence of an unseen third body in the sys- 
tem, orbiting at a great distance the closer central pair. A 
discussion of this hypothesis is beyond the scopes of the 
present paper. We may however point out that in our spec- 
tra when the Balmer lines show a clear split and the S/N 
ratio is high, a weak absorption is visible. The Ha com- 
ponent does not correspond to known telluric absorptions 
in this region, and the coincidence between the various 
Balmer lines favour an interpretation as a real stellar line. 
To the benefit of future Investigators of SS Lac we report 
in Table 3 the heliocentric velocity of this weak component 
as measured on our spectra (same weighting as described 
in Sect. 2), that could in some way be related to the unseen 
third body. 
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